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Abstract 

We have investigated constraints on the coupling between dark matter and the interacting 
Chaplygin gas. Our results indicate that the coupling constant c between these two entities can 
take arbitrary values, which can be either positive or negative, thus giving arbitrary freedom to 
the inter-conversion between Chaplygin gas and dark matter. Thus, our results indicate that the 
restriction < c < 1 on the coupling constant occurs as a very special case. Our analysis also 
fN| supports the existence of phantom energy under certain conditions on the coupling constant. 
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1 Introduction 
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It is well known that the expansion of the universe is accelerated. It has been confirmed by 
numerous observations taken by various scientific groups across the globe using WMAP 
[Tj, distant supernova type la data [2j [3], large-scale structure and galaxy distribution p|] 
and the gravitational lensing phenomenon of high-redshift galaxies [5j . These observations 
clearly suggest that the universe is spatially flat and is dominated by some sort of vaccum 
energy having negative pressure commonly called 'dark energy'. This energy has been 
interpreted in various forms like the cosmological constant jSJ [7], quintessence models 
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based on the ideas of a spatially homogeneous and time-dependent scalar field [E] , phantom 
energy [91 [101 HU [391 H2]) quintom models [HI [151 [16], k-essence [17], holographic dark 
energy [18] and the Chaplygin gas (CG) [19], [20] (see also |21j for a recent review on dark 
energy). 

The CG is represented by an equation of state (EoS) of the form p = —A/p, where 
A is a constant parameter |22j. The CG gives rise to a simple cosmological model that 
interpolates between the earlier matter- (or dust-) dominated to the later dark energy 
dominated phase of the universe. Due to its effectiveness in explaining the evolution 
of the universe, several generalizations of CG have been proposed in the literature [231 
[2U [251 12S] . The observational evidence in support of cosmological models based on the 
CG-EoS is also very encouraging [271 EH] • The CG also possesses the property of giving 
accelerated expansion even if it gets coupled with other scalar fields like quintessence or 
dissipative matter fields [20] • It also yields traversable wormhole solutions to the Einstein 
field equations if the pressure and density of CG violates the null energy condition [30] . 
Besides its various useful implications in cosmology, CG has the drawback of producing 
oscillations or exponential blow up of dark matter power spectrum which is inconsistent 
with observations [21] ■ Similar results are obtained in later generalizations of CG |32j . 
But later it was proved that such oscillations can be avoided and structure formation can 
proceed (which is strongly supported by dark matter) if the phantom-like dark energy is 
excluded, thereby proceeding with only dark matter and dark energy [33] ■ It was further 
suggested that CG behaved like a passive background in the early evolution of the universe 
and that only dark matter leads to nonlinear growth of structures but later the evolution 
is dominated by CG [31]. The inhomogeneities generated by dark matter were stabilized 
by the CG which is compatible with the observations [35]. Later studies on supernovae 
data put constraints on CG leading to cosmological models based on CG to behave just 
like the cosmological constant [36J. 

Modern cosmology is plagued with numerous theoretical and observational problems: 
among them is the cosmic- coincidence problem which can be stated thus |37j : why are 
the energy densities of matter and dark energy almost of the same order at present? In 
the standard cosmological model, the ratio of the energy densities of matter and dark 
energy should fall rapidly as the universe expands, but observationally the corresponding 
ratio turns out to be almost constant or minutely fluctuating around unity, a phenomenon 
commonly called the 'soft coincidence'. It leads to the possibility that energy might be 
exchanged to keep such a delicate balance in the densities. This interaction is generally 
studied in the models so-called 'interacting dark energy' [381 ESI BH HI]- Unified models 
based on dark matter and CG have been widely investigated (see [33] and references 
therein) but the fundamental question dealing with the interaction between these two 
entities is not satisfactorily answered and requires further investigation. A cosmological 
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model based on the interacting CG had been proposed [32] to investigate this interaction. 
This model yields the result that the universe is to cross the phantom divide i.e. the 
transition from the state w > -1 to w < -1 or more simply u — — 1, which is not 
possible in the models based on pure CG. Furthermore, this leads to the scaling solutions 
of the cosmological dynamical system, which helps in explaining the coincidence problem 
effectively. Moreover, it is used in cosmological models to investigate dissipative effects 
of van der Waal's fluid and dark energy [33] • In fact it is later suggested that without 
interaction with other species, any hydrodynamical or k-essence like model in general 
relativity cannot cross uj = — 1 [33]- The interacting CG also yields stable scaling solutions 
of Friedmann-Lemaitre- Robertson- Walker (FLRW) equations at late times of the universe. 
This model was later extended to the case of an interacting generalized Chaplygin gas 
(GCG) [32]- There are some proposals that this interaction can be observed if cubic 
corrections are provided to the Hubble law, measured by distant supernovae of type la 
[3S]- It is worthwhile to understand the role of the coupling constant in the interacting 
models which we have investigated using a modified Chaplygin gas (MCG). 

The outline of this paper is as follows: In the next section, we model our dynamical 
system on the pattern of [32] and determine the critical points corresponding to that 
system. In the third section, we perform stability analysis corresponding to each critical 
point. In the fourth section, we perform analysis to determine constraints on the coupling 
constant. Finally we present conclusion of our paper. 



2 Modeling of dynamical system 

We start by assuming the background to be a spatially homogeneous and isotropic FLRW 
spacetime, 



ds 2 = -dt 2 + a 2 (t) 



dr 2 



+ r 2 {d6 2 + sin 2 



1 — kr 2 

where a(t) is the scale factor and the curvature parameter k = —1, 0, +1 describes spatially 
open, flat or closed spacetimes. It is assumed that the spacetime is filled with the a 
two-component fluid, namely dark matter and dark energy. The corresponding energy- 
momentum tensors are specified by 

Tap ^ Pdrn^ u}^vl T^lu ' (Pmcg Vmcg)^ fx^p Vmcg^ /i^u • (2) 

Here u'^ and is the comoving four-velocity of dark matter and dark energy respectively. 
The notations dm and meg corresponds to dark matter the dark energy, respectively. In 
our model, the dark energy is specified by the modified Chaplygin gas EoS [47J 

B 

Pmcg -^-Pmcg ~ ; \<J) 

Pmcg 
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where A and B are constant parameters and < a < 1. The MCG reduces to GCG if 
A = and to the CG if furthermore a = 1. It reduces to the standard linear EoS for a 
perfect fluid if B — 0. Recently it has been deduced, using latest supernova data, that 
models with a > 1 are also possible [H]. For our analysis the positivity of a is sufficient. 
In our further discussion, the MCG and dark energy are used interchangeably. 

The density evolution of MCG is given by 

( B C 

Pmcg ~\ l + A + a 3(l+A)(l+a) ) ' W 

where C is the constant of integration. Note that Eq. (4) holds only when the interaction 
is absent. The equations of motion corresponding to FLRW spacetime filled with a two- 
component fluid are 

K 2 

H = ~^r(Pmcg 4~ Pmcg 4~ Pdm): (5) 

K 2 

H 2 = —{Pmcg + Pdm)- (6) 

Here k 2 = 8ttG is the Einstein gravitational constant and H = H(t) is the Hubble 
parameter. We assume k = 0, representing a flat model of the universe. Furthermore, the 
energy conservation for the two-component perfect fluid is obtained from 

= V U (T^ + T^) = 0. (7) 

Here refers to the covariant derivative with respect to x u coordinate. Eq. (7) yields 

Pmcg 4" pdm + 3H(p mcg + p mcg + p dm ) = 0. (8) 

Due to interaction, the energy will not independently be conserved for the interacting 
components, and therefore 

V u T imcg ) = _ Q ^ VUT (dm )=Q ^ (g) 

Here is the interaction term that corresponds to energy exchange between dark energy 
and dark matter. Solving Eqs. (9) using (2), we obtain the so-called energy-balance 
equations corresponding to MCG and dark matter as: 

Pmcg 3i? (Pmcg Pmcg) -Q, (io) 
p dm + 3Hp dm = Q. (11) 

The function Q = Q t , p = t has dependencies on the energy densities and the Hubble 
parameter, i.e. Q(Hpd m ), Q(Hp mcg ) or Q(Hpd m , Hp mcg ) jl9]. Because of the unknown 
nature of both dark energy and dark matter, it is not possible to derive Q from first 
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principles. In order to deduce a reasonable Q, we may expand like Q(Hp dm , Hp mcg ) ~ 
a dm H p dm + a mcg H p mcg . Since the coupling strength is also not known, we may adopt just 
one parameter for our convenience; hence we take a dm = a mcg = c [50]. We here choose 
the following coupling function Q given by [39] : 

Q = 3Hc(p mcg + p dm ). (12) 

The choice of Q is completely arbitrary but care must be taken that it must satisfy 
the energy conservation (see Ref. [52] for various other forms of Q). Here c is the 
corresponding coupling constant (also called the 'transfer strength') for the interaction. 
To study the dynamics of our system, we proceed by setting 

x = lna = — ln(l + z), (13) 

where z is the redshift parameter. Moreover, the density and pressure of MCG can be 
expressed by dimensionless parameters u and v as follows: 

2 2 

r^. Pmcq rJ Pmcq K Pmcq /t a\ 

tt=n ^ = ^r = w v = ^ip- (14) 



The EoS parameter uj is conventionally defined as 

Pmcg 
Pmcg 



U(X) EE (15) 



which becomes 

u(x) = -. (16) 
u 

The density parameters of MCG and dark matter are related as 

ft} 

^dm ^ ^mcg 1 U. (17) 

Since Q dm ~ 0.3 [2J, it constrains u G (0, 1) for a flat universe. The following system of 
differential equations governing the dynamics of our cosmological model is determined by 
using the above equations: 

du , , 

— = -3c-3v + 3uv, 18 
ax 

— = -3 -+(A--)(l + a) (u + v + c) + 3v(l + v). 19 
ax lu V u J 



Note that for A = 0, the above system reduces to the case for interacting generalized 
Chaplygin gas |15] . By equating Eqs. (18) and (19) to zero, we obtain the three critical 
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points (u ic , Vi C ), with i — 1, 2, 3, given by 

u lc = 1-c, (20) 

uic = -1, (21) 
1 VA + 4c 

U2c = 2-"W' (22) 

w 2c = ^(A-y/Ay/A + Ac), (23) 

1 / yJA + 4c\ , , 

«3 C = „ 1 + » (24) 



2 V VA 

v 3c = ~{A + VaVA + Ac). (25) 

For the critical points to be real valued, we require A + 4c > 0. Notice that the first 
critical point is the same as discussed in [41j] . There it was proposed that the coupling 
constant c G [0,1]. Since there is a transition from CG to dark matter (i.e. c — > 1) as 
the universe evolves, it implies that the future universe will contain only dark matter and 
might have no trace of CG. Our analysis in the next two sections suggests that c cannot 
necessarily be restricted in the range < c < 1 and can take values outside this range. 



3 Stability analysis 

To perform a stability analysis of our dynamical system, we linearize the system of equa- 
tions (18) and (19) about the critical points to get 

dSu 

— - = 3v c 5u + 3(-l + u c )5v, (26) 
ax 

dSv 3 

— - = -[av c {c + v c ) + Au 2 c {l + a)]5u (27) 

3 

H — [(c + 2v c )a + u c (l + 2v c + a- A(l + a))]5v., 

The eigenvalues of the above dynamical system (26) and (27) corresponding to the three 
critical points Eqs. (20 - 25) are 
3 

Ai = — — (2 + A-2c- Ac + a + Aa- Aca (28) 

2{c — I) 

+ + A(c - l))(c - 1) 2 (1 + a) + (2 - 2c + a - A(c - 1)(1 + a) 2 )), 

_3 

H X = — tt(-2- A + 2c + Ac- a- Aa + Aca (29) 

L \ C 1 ) 



+V4(-1 + A(c - l))(c - 1) 2 (1 + a) + (2 - 2c + a - A{c - 1)(1 + a) 2 )), 
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A 2 = -^(-2 - 2a - A(l + a) - VAVA + 4c(3 + a) + V2[{2A 5/2 a + 2(1 + a) 2 

+2\fA^A + 4c(l + a) 2 + 2A(l + a)(-l + c - 2vM. + 4c + a + ca) (30) 
+A 2 (1 - 2\J A + 4ca + a 2 ) + A 3/2 (4 - y/A + Ac + 4a - 8ca + vM + 4ca 2 ))] 1/2 ), 

\H = ^(2 + 2a + A(l + a) + V^vM. + 4c(3 + a) + v / 2[(2A 5/2 a + 2(1 + a) 2 

+2\fA^A + 4c(l + a) 2 + 2A(1 + a)(-l + c - 2vM. + 4c + a + ca) (31) 
+A 2 (1 - 2V A + 4ca + a 2 ) + A 3/2 (A - y/A + 4c + 4a - 8ca + + 4ca 2 ))] 1/2 ), 

A 3 = (-2 - 2a - A(l + a) + V^vM + 4c(3 + a) + v / 2[(2A 5/2 a + 2(1 + a) 2 

-2v / Av / A + 4c(l + a) 2 + 2A(1 + a)(-l + c - 2V A + 4c + a + ca) (32) 
+A 2 (1 + 2vM + 4ca + a 2 ) + A 3/2 (4 + y/A + 4c + 4a + 8ca - + 4ca 2 ))] 1/2 ), 



/i 3 = |(2 + 2a + A(l + a) - V^vM. + 4c(3 + a) + v / 2[(2A 5/2 a + 2(1 + a) 2 

-2y/Ay/A + 4c(l + a) 2 + 2A(1 + a)(-l + c - 2vM + 4c + a + ca) (33) 
+A 2 (1 + 2y/A + Aca + a 2 ) + A 3/2 (4 + y/A + 4c + 4a + 8ca - + 4ca 2 ))] 1/2 ). 

It is easy to show that the real parts of the eigenvalues (Ai, /ii) are negative, while for 
(A 2 , /i 2 ) and (A 3 , /i 3 ), the real parts are all positive. Hence a stable stationary attractor 
solution is possible through the first critical point only. It can alternatively be proved 
using the deceleration parameter q. Note that the acceleration in the late evolution of the 
universe arises when 

5 =-Ji<-l, (34) 

which holds for the first critical point (ui c , Vi c ) only, since qi = —1. While for (u 2c , v 2c ) 
and (u 3c ,v 3c ), we require 



1 

Q2,3 - ^ 



1 + ^(AtVaVA + 4c) 



< -1. (35) 



As A + 4c > 0, the inequalities in Eqs. (35) do not hold and hence the accelerated- 
expansion solution is not obtained from the second and third critical points. Hence the 
valid attractor solution is obtained from the first critical point only. We shall, henceforth, 
deal with the first critical point only. 

As shown in figure 1, the first critical point (ui c , Vi c ) is the stationary attractor solution 
for the interacting modified Chaplygin gas with the coupling constant fixed at c = 0.5. 
The parameter A can assume values in the range —0.35 < A < 0.025 |54l : while we 
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choose A = 0.025 for our numerical work. Also note that if the parameter a < 0, 
then it yields a polytropic equation of state for dark energy but for the MCG, we take 
a = 0.004. It is evident that all the solutions of the dynamical system with four different 
initial conditions converge to the same final state. As q\ = —1, the first critical point 
gives rise to an accelerated-expansion solution of the universe which is consistent with the 
observations. 

Moreover, the attractor solution corresponding to (ui c , V\ c ) is also possible if c takes 
values outside the usual considered range of < c < 1. In figures 2 and 3, the parameter 
c is given values 1.7 and —1.5, respectively, with the same initial conditions. Curiously, 
all the four solutions converge to the same single final state. It draws to the fact, that at 
least theoretically, the coupling constant c can take values outside the interval [0,1]. This 
result is further deduced in the next section using a different formalism. 



4 Constraints on coupling constant 



We can determine the constraints on the coupling constant c by using the first critical 
point of our dynamical system. For this purpose, we shall adopt the formalism of Guo 
and Zhang [35] . We define new parameters corresponding to MCG and dark matter by 



Pmcg Pmcg 



and 



Idr, 



Pmcg 

Pdm + Pd in 
Pdm 



(36) 



(37) 



Note that jdm = 1 since pd m = 0. Moreover, the parameter 7 mc9 will be determined 
corresponding to the first critical point. To find how the density ratio R evolves with 
time, we differentiate it with respect to t to get 



R 



dR 
~dt 



Pmcg 



Pdn 
Pdn 



Pmcg 
Pmcg 



(3? 



Using Eqs. (10) and (11), Eq. (38) becomes 

R = R 



1 h 3# (7 mc5 



Pdm Pmcg 

Using Eq. (12) in (39), we get after simplification 

R = 3H[c(l + Rf + R{% 



mi g 



!)]■ 



(39) 



(40) 
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In order to get stationary solutions, we solve for R = to get 



Rf= 1 2 ^ mcg - 1± y( 1 2c mcg ~ 1 ) _L (41) 

Now, to get real valued solutions, we require ( 1 ~^ lcg — l) 2 — 1 > 0, which yields 

(i^-2)(l^=)>0. (42) 

The above inequality holds if the quantities in the brackets are either both positive or 
both negative. We shall take c to be a free parameter which can take values other than 
zero. 

Case (1) 

Assume both quantities in the brackets in (42) to be positive, i.e. 

1 - 7mcg - 2 > 0, 1 - lmcg > 0. (43) 
2c ~ ' 2c y J 

Case (la) 

Now take c > 0; thus Eq. (43) gives 

lmcg — 1 4c, lmcg — 1; (^4) 

which yields 

lmcg < 1 - 4C. (45) 

Case (lb) 

If c < 0, then Eq. (43) yields 

lmcg > 1 - 4C, l mcg > 1, (46) 

which implies 

lmcg > 1 - 4C. (47) 

Case (2) 

Now take both quantities in the brackets in Eq. (42) to be negative, i.e. 

1 ~ lmcg - 2 < 0, 1 - lmcg < 0. (48) 
2c ~ 2c K J 

Case (2a) 
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Take c > 0; after solving Eq. (48), which gives 



which yield 

Case (2b) 

If c < 0, then Eq. (48) yields 



Imcg — 1 4c, ^mcg — lj 
Imcg — 1 • 

Imcg ^ 1 4c, 'Jmcg — lj 



(49) 
(50) 

(51) 



which results in 

7 mcs < 1. (52) 

Now we shall use the definition 7 mcs = 1 + ui = 1 + v\ c j u\ c in each of the above four 
cases. 

Case (la) 

Using uj\ = ^ = ^ in Eq. (45), we have (1 — 2c) 2 > ,which is satisfied for all 
values of c. Notice that from Eq. (44), we have an additional constraint c < 1; therefore, 
< c < 1, which is the range usually considered for c in the literature. 

Case (lb) 

Using uoi in Eq. (47) we get (1 — 2c) 2 < 0, which is satisfied only for c = 1/2. Since 
c < 0, we do not have an acceptable solution. 

Case (2a) 

Here for cui, Eq. (50) implies ^ > 0, which holds for all c > 1. Apparently it implies 
that the coupling constant between MCG and dark matter can take arbitrary value; thus, 
the mutual interaction can be more dynamic. It yields arbitrary freedom for the conversion 
of MCG into dark matter. 

Case (2b) 

For u>i, Eq. (52) implies < 0, which is viable if c < 0. Thus, the coupling constant 
can take arbitrary negative values. This case apparently supports the conversion of dark 
matter into MCG with arbitrary coupling. Note that iO\ represents the EoS of phantom 
energy (wj < -1). 



10 



5 Conclusion 



We have investigated the possible interaction between dark matter and the Chaplygin gas 
and we deduced that the coupling constant involved can take values outside the range 
usually considered, < c < 1. This range arises as a special case in the Case (la). Our 
analysis suggests that c can take arbitrary positive or negative values. If c > 1, as in 
Case (2a), then it supports the conversion of MCG into dark matter. Conversely, if c < 
as in Case (2b), it allows for the conversion of dark matter into MCG. It also supports 
the existence of phantom energy through the final case. Moreover, the present work may 
serve as the generalization of the earlier work by Zhang and Zhu [12] for the interacting 
Chaplygin gas and by Wu and Yu [45J for the interacting generalized Chaplygin gas. 

In a recent investigation, Feng et al [33] have presented observational constraints on the 
coupling parameter and have deduced that small positive values for c are most probable. 
This conclusion is drawn in order to alleviate the cosmic-coincidence problem. Also the 
negative values of c are excluded to avoid the violation of the second law of thermodynam- 
ics [22] • We have deduced from our analysis that c has no such theoretical constraints, 
and the usual choice [0,1] is not a true range for the coupling parameter. 
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Figure 1: The phase diagram of the interacting modified Chaplygin gas model with c = 0.5. 
The model parameters are fixed as A = 0.025 and a = 0.004. The curved lines from left to right 
correspond to the initial conditions u(— 2) = 1.2, v(— 2) = —0.2 (green); u(— 2) = 1.3, v(— 2) = 
-0.3 (blue); u{-2) = lA,v{-2) = -0.4 (red); u{-2) = 1.5,v(-2) = -0.5 (black). 
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Figure 2: The phase diagram of the interacting modified Chaplygin gas model with c = 1.7. 
The model parameters are fixed as in Fig. 1. The curved lines from left to right correspond to 
the initial conditions as given in Fig. 1. 
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Figure 3: The phase diagram of the interacting modified Chaplygin gas model with c = —1.5. 
The model parameters are fixed as in Fig. 1. The curved lines from left to right correspond to 
the initial conditions as given in Fig. 1. 
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